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We present updated preliminary results for the nucleon electromagnetic form factors for non- 
perturbatively &{a) improved Wilson fermions in Nf = 2 QCD measured on the CLS ensembles. 
The use of the summed operator insertion method allows us to suppress the influence of excited 
states in our measurements. A study of the effect that excited state contaminations have on the 
i Q 2 dependence of the extracted nucleon form factors may then be made through comparisons of 

the summation method to standard plateau fits, as well as to excited state fits. 
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a[fm] 


L/a 


L[fm] 


m^[MeV] 


no. meas 


5.2 


0.079 


32 


2.5 


312, 363, 


1696, 796, 










473, 603 


1060, 576 


5.3 


0.063 


32 


2.0 


451, 606, 649, 


1344, 1296, 278, 






48 


3.0 


277, 324 


1000, 796 


5.5 


0.050 


48 


2.4 


430, 536 


600, 600 



Table 1: Details of the lattice ensembles used in this study, showing j3-values, lattice spacing a, lattice 
extent L (where T — 2L), pion mass m K and total number of measurements. 



1. Introduction 

Baryon form factors are central observables of hadronic physics and provide details of the 
object's distribution of charge and magnetisation. Currently, lattice simulations fall short of the 
accuracy achieved by experiment. Furthermore, the simulations of the nucleon electromagnetic 
form factor fail to reproduce experimental results [1,2]. It is therefore important to ensure that sys- 
tematic effects are under control in lattice simulations. To study one such effect, we employ three 
separate methods (described in section 3) that account for excited state contributions with varying 
degrees of rigour to check the control of systematic errors in the extraction of the form factors. This 
work provides an update on results previously presented in [3] and follows the methodology used 
in a recent study of the nucleon's axial form factor [4]. Similar methods have been used in [5]. Our 
simulations use non-perturbatively &{a) improved Wilson fermions in Nf = 2 QCD, measured on 
the CLS ensembles. Table 1 provides details of the lattice ensembles. 

The matrix element of a nucleon interacting with an electromagnetic current V 1 = \\f(x)y^ y(x) 
may be decomposed into the Dirac and Pauli form factors F\ and F2: 



(N(p , ,s')\V^\N(p,s)}=u(p',s') 



ZrriN 



u(p,s), (1.1) 



where u(p,s) is a Dirac spinor with spin s, and momentum p, 7^ is the Dirac matrix, and a^ v = 
hlYft-iYv]- Also Q 2 = —{Epi — E p ) 2 + q 2 where q = p' — p. These form factors are related to the 
Sachs form factors, Ge and Gm, 

<**>-**,-***,. <«*>-**>+**,. 0, 

that are measured in scattering experiments via the differential cross section described by the 
Rosenbluth formula. The form factors may be Taylor expanded in the momentum transfer Q 2 

Gx(Q 2 ) = Gx(0) (l + \{r 2 )Q 2 + ^(g 4 )) , (1.3) 
from which the charge radii of the nucleon may be determined: 

6 dG x (Q 2 



(4) 



Gx(Q 2 ) dQ 2 



(1.4) 

2=0 



where X = E, M. For the conserved current, Ge(0) = 1 and Gm(0) = pi, where pL measures the 
magnetic moment in nuclear magneton units e/(2m#). 
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2. Lattice formulation 

The calculation of the form factors requires a ratio of correlation functions, for which we use 
(for the case p' = 0) 



R r Mt,t s ) = WJ^-WWt), (2 .d 

r " C 2 (0,t s ) V C 2 (0,t s -t)C 2 (q,t)C 2 (q,t s y 

C 2 (P,t) = Y,(r a , a J a (x)J a ,(0)}e- r P rx , (2.2) 

X 

C 3 ^{q,t,t s ) =Y i {T a , a J a {x,t s )ff yil {y,t)J al {0))e- rq -y. (2.3) 

x,y 

This ratio was found to be the most effective in [6]. C 2 (p,t) and Cs.y^(q,t,t s ) are two- and three- 
point functions respectively. J a (x) is a suitably chosen interpolating operator with the correct 
quantum numbers to create a nucleon, and T aa i is a projection matrix used to give the interpolating 
fields the correct parity. We choose to polarise the nucleon in the z-direction, for which T = 
2~(1 + 7o)(l +'7573)- This choice allows us to extract both Ge and Gm, whereas for an unpolarised 
nucleon, F = g(l + %), one may only determine Ge- Further to this, for improved statistics we 
average over polarisation in the positive and negative z-direction as well as over the nucleon and 
anti-nucleon state. 

Due to the insertion of the operator at time t, the three point function is a more computationally 
demanding object. We use the 'fixed sink method' for its calculation, which fixes the final and 
initial states, but allows both the operator and momentum transfer to be chosen without the need 
for additional inversions [7]. For this study, we consider both the local and conserved vector current 
for the operator . The latter is defined as 



scon 
i 1 



1 



(*) = ~ ( V(* + a A ) ( 1 + 7m ) u l M ¥(x) - V(x) ( 1 - 7 M ) Up (x) \j/(x + a fx ) ) (2.4) 



where y = u,d. In principle, we are able to determine the form factors for the proton, neutron, iso- 
scalar and iso-vector depending on the linear combination of contributions from the quark correla- 
tion functions. We concentrate on the iso-vector, which has the advantage that quark-disconnected 
diagrams cancel. 

To improve the overlap of the interpolating operators with the nucleon, we used Gaussian 
smearing at both the source and sink supplemented by HYP smeared links. Whilst this is a great 
improvement over point sources, the extraction of plateaus at non-zero momenta can still prove 
difficult. A poster presented at this conference addresses this problem for two-point functions, 
whereby using a generalisation of Gaussian smearing to form an anisotropic wavefunction results 
in a reduction of the noise-to-signal ratio at non-zero momenta for the pion and nucleon [8, 9]. 
We are currently extending this study to include three-point correlation functions with the aim of 
applying the technique to the extraction of the vector form factors. 

Ge and Gm may be extracted from eq. (2.1) at large time arguments: 



\/ A1 ^Ge(Q 2 ), R r M = EijPjJ 2e{ e + M ) Gm{Q2) 1 = 1)2 - (2 - 5) 
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Figure 1: Left panel: Plateau fits to Ge for several t s for the smallest non-zero momentum transfer, Q 2 . 
Excited state contamination effects for the smallest t s are clearly visible (the asymmetric distribution of 
points is a remnant of the non-zero momentum transfer). Right panel: Summation method for Ge- Both 
panels are shown for our lightest ensemble (m K — 277 MeV). 



3. Systematics of extraction 

Correlation functions must have reached their asymptotic behaviour for a reliable and unbiased 
determination of the form factors. However, we observe exponentially decaying excited states 
from both the source and sink. Therefore, simple plateau fits (left panel fig. 1) show a trend of 
higher values for small source-sink separations, i.e. for decreasing t s . To control systematics, it is 
important to take these excited states into account. Contributions to the ratio from the ground and 
excited states may be factorised 

R(q,t,t s )=R o (q,t,t s )(l + 0{e- At )+0{e- A '^)) (3.1) 

where A and A' are the energy gaps of the initial and final nucleons respectively. With the assump- 
tion A = A' = 2m K l we take the excited states into account using a fit to 

R{q,t,t s ) =GE,M + bie- At +b 2 e- Ait ^+b3e- At \ (3.2) 

An alternative without the need for the assumption A = A' = 2m n uses summed operator insertions 
[10]: 

S(t s ) = f,R(q,t,t s ) -> c(A,A') +t s (G EM + 0(e- At <) + . (3.3) 

This allows the form factors to be extracted from the slope, from computing S(t s ) for several t s 
(right panel fig. 1). The excited states should be more suppressed for this method than for a fit to 
eq. (3.2), because t s > t, (t s — t). 

'This is not strictly true when there is a momentum transfer, however we find this to be a small effect as the data 
is well described by eq. (3.2). Further to this, the results agree with the summation method, which does not require an 
assumption for the energy gap. 
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Figure 2: Results for our most chiral (m K = 277 MeV) lattice. The top two panels show the Q 2 dependence 
for Ge and Gm- The bottom panel shows a determination of the magnetic moment, ji. This is shown for the 
three extraction methods. 

4. Form factor Q 2 dependence and chiral behaviour 

For the remaining discussion we concentrate on the conserved current as this removes the need 
for any renormalisation of the lattice operators; however, we note that a comparison between the 
local and conserved current provides a check of the renormalisation factor, which we find to be in 
agreement with other work (e.g. [11]). In order to model the Q 2 dependence of the form factors 
(shown in fig. 2 for our most chiral analysed ensemble) a dipole ansatz is adopted 

GWO) 



Gem(Q 



M l.M 



(4.1) 



from which the radii is extracted. We may also obtain the magnetic moment pL from Gm{0) and 
also from the ratio 

G M {Q 2 ) 



lim 



(4.2) 



r e^oG £ (2 2 )' 

based upon the phenomenological observation that the form factors Ge and Gm have very similar 
radii. For our determination, the two methods are in good agreement (fig. 2). 

This calculation was performed on the ensembles listed in table 1 to check both finite volume 
and discretisation effects. The lattice ensembles cover a range of pion masses (277 to 649 MeV). 
In order to perform an extrapolation in the pion mass to the physical point, we model the chiral 
behaviour of the radii using a HB^PT inspired fit [12], 

{rj) = c l + c 2 log(ml). (4.3) 
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We are currently investigating the effect of other fit forms on the result, including linear fits and 
covariant B^PT fits as well as applying cuts to the pion mass range [13, 14, 15, 16], so as to have 
a comprehensive picture of the systematic effects. We also apply this to (r 2 ) and fc = /I — 1. 

Thus far we have sketched the methodology commonly employed to extract form factors on 
the lattice. Using this, we study the systematics of the extraction by separately employing the 
three analyses described above: a plateau fit (for the largest t s ), a simultaneous excited state fit to 
t and t s and the summation method. A comparison of the three methods in fig. 2 shows that whilst 
all three methods agree within statistical errors, a systematic trend for better agreement between 
the two methods that account for excited state effects is observed. This is echoed in all of our 
ensembles. An earlier study by our group on gA [4] demonstrates a clearer systematic dependence. 
However, we note that gA is a somewhat cleaner quantity to extract on the lattice (i.e. Q 2 = 0). 
The electromagnetic form factor analysis requires an improvement in the statistical errors and the 
addition of smaller pion masses to further check the systematic behaviour and as such the results 
in this proceeding should be considered as preliminary. Fig. 3 compares the determination of the 
Dirac radius (rj) from an excited state fit with the plateau fit (the summation method agrees well 
with the excited state fit but with increased statistical errors). We observe an increase in the central 
value for all but one ensemble and also that the difference increases as the pion mass decreases (as 
one would expect, excited states should contribute more for more chiral m K ). This results in a larger 
radius that is closer to the experimental result. If we restrict the fit to the four most chiral ensembles 
and perform a straight line fit, a result compatible with experiment is achieved for the excited state 
method for (r 2 ). We have also looked at fc and (r 2 ), for which we see a similar behaviour. These 
are however all at the preliminary stage. 

5. Conclusions and outlook 

We have presented preliminary results for the nucleon vector form factors with a particular 
emphasis on the potential for systematic errors from excited state contaminations. We account 
for this through a comparison of three extraction techniques which account for excited states with 
varying degrees of rigour. Whilst we are, so far, unable to see any effects outside of errors between 
the three methods, we do see an improved agreement between the two methods that account for 
excited state effects in both the Q 2 behaviour and in the chiral extrapolations of the form factors, 
and consider this to be an indication of a systematic effect. Furthermore, there is a trend for a 
steeper gradient to be extracted from these results and thus such a systematic could lead to better 
agreement between lattice results and experiment. The results indicate the importance of excited 
states and that they should be a consideration in studies of potential systematic effects. 

The large statistical errors for all methods highlight the need for greater statistics and for more 
chiral points to be added, where excited states should, in principle, have a greater effect. Monte 
Carlo ensembles exist for more chiral points (the lightest available is approximately 200 MeV), but 
are yet to be analysed. 
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Figure 3: Chiral extrapolation of the Dirac radius to the physical point (vertical black line) using {r\ ) = c\ + 
C2log(m^) using the largest t s plateau-fit (blue) and excited state fit (red) for the extraction. The horizontal 
black line shows the experimental result. The different symbols indicate the lattice spacing (see legend) and 
the leftmost two points show our extrapolated values. 
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